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Non-trivial critical models in 2D with central charge c = are described by Logarithmic Con- 
formal Field Theories (LCFTs), and exhibit in particular mixing of the stress-energy tensor with a 
"logarithmic" partner under a conformal transformation. This mixing is quantified by a parameter 
(usually denoted b), introduced in [V. Gurarie, Nucl. Phys. B 546, 765 (1999)]. The value of b has 
been determined over the last few years for the boundary versions of these models: 6 perco = — | 
for percolation and b po iy = § for dilute polymers. Meanwhile, the existence and value of b for the 
bulk theory has remained an open problem. Using lattice regularization techniques we provide here 
an "experimental study" of this question. We show that, while the chiral stress tensor has indeed 
a single logarithmic partner in the chiral sector of the theory, the value of b is not the expected 
one: instead, b — —5 for both theories. We suggest a theoretical explanation of this result using 
operator product expansions and Coulomb gas arguments, and discuss the physical consequences on 
correlation functions. Our results imply that the relation between bulk LCFTs of physical interest 
and their boundary counterparts is considerably more involved than in the non- logarithmic case. 

PACS numbers: 11.25.Hf, 02.30.1k, 64.60.ah 

Introduction. Two-dimensional (2D) logarithmic con- 
formal field theories (LCFTs) with central charge c = 
describe the scaling limit of geometrical problems such 
as polymers (self-avoiding walks) or percolation (see 
e.g. [H). They also play a fundamental role in the study 
of phase transitions in systems with quenched disorder 
|2j-|4|: for instance, the long sought-after conformal field 
theory describing the plateaux transition in the integer 
quantum Hall effect should be of this type. Other phys- 
ical applications range from non-equilibrium systems |6[ 
to aspects of the AdS/CFT correspondence (see e.g. [3]) 
and description of super-symmetric sigma models beyond 
the topological sector Q . 

Progress in this field since the pioneering papers 0, H3] 
has been difficult. Chiral (i.e., boundary) LCFTs are 
slowly getting under control thanks to recent progress 
on the abstract study of non semi-simple modules of the 
Virasoro algebra 



ries at c = (see e.g. [17!, Il8j ) often exhibit unwanted 



-the so-called staggered modules [11 [ . 
Consistent results have also been obtained from lattice 
approaches 
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131 ] . It is now well-accepted that chiral 
LCFTs are characterized by a complicated structure of 
such staggered modules, and an infinity of numbers called 
indecomposability parameters, or logarithmic couplings. 
One of these numbers, called b [3, H[, is of particular 
interest as it encodes the structure of the module of the 
stress-energy tensor in c = theories. Numerical 141 115| 
and theoretical fl(| arguments lead to the values & per co = 
— I for percolation [Q —> 1 state Potts model] and 6 po iy = 
I for dilute polymers [n — > dilute O(n) model]. 

While chiral LCFTs are thus fairly well understood, 
much remains to be done to understand the structure 
of bulk theories. Attempts to construct non-chiral theo- 



features such as degenerate or non SL(2,C) invariant 
ground states which should not occur in, for instance, 
percolation. 

Part of the progress in the chiral case originates from 
considering lattice models with open boundary condi- 
tions, and observing that the indecomposable features 
of the chiral (Virasoro) algebra appearing in the scaling 
limit are similar (for more accurate, mathematical state- 
ments, see (lij ) to those occurring, in finite size, on the 
lattice. This suggests that eventually, LCFTs may be 
solved by a careful exercise in the representation theory 
of the (associative) algebras satisfied by the local energy 
terms (such as the Temperley-Lieb algebra [2(3|). The 
non-chiral or bulk case corresponds, on the lattice, to 
periodic boundary conditions, which is rather difficult 
mathematically [2l| . One can nevertheless use lattice al- 
gebraic techniques to investigate aspects of the simplest 
modules under the full left and right Virasoro algebras 
present in this case. The full results will be discussed 
elsewhere [22| , but, as a first step, we present in this Let- 
ter the results of the corresponding numerical measure 
of the parameters b for both percolation and polymers 
with periodic boundary conditions (PBC). The outcome 
is quite unexpected: both theories turn out to have the 
same 6, which moreover differs from that of the open 
(chiral) case. 

b parameter. We begin by defining the b number, first 
introduced by Gurarie [3|-|5|. We consider the (non- 
chiral) LCFTs that describe the scaling limit of polymers 
or percolation. Within these theories, the stress-energy 
tensor T{z) = L-2I is mixed into a Jordan cell with its 
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FIG. 1: Example of dense (left panel) and dilute (right panel) 
loop configurations. We use periodic boundary conditions in 
both directions so that the loops actually live on a torus. We 
also show the corresponding FK clusters for the dense loop 
model. Each closed loop has a weight n. 



logarithmic partner t(z,z) such that in the basis (T, t), 
the generator of the scale transformation reads 



Lr 



u = 



2 2 
2 



(1) 



We will show later (see also [23|) that the non- 
diagonalizablc term can be decomposed as Lgt = 2t + T 
and Lgt = T, so that, unlike T(z) which is purely holo- 
morphic, the field t(z, z) has a non-trivial antiholomor- 
phic part dt ^ 0. Using this normalization, we define the 
number b as 



(T\t), 



(2) 



where we have used the bilinear Virasoro form. Note also 
that Lit — and L2t = bl. Using these relations, global 
conformal invariance fixes Q the form of the two-point 
functions to be 



(T(z)T(O)) 
(T(z)t(0,0)) 

(t(M)*(0,0)) 
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(3a) 
(3b) 

(3c) 



where 9 is an irrelevant constant that can absorbed into 
a redefinition of the fields, whereas 6 is a fundamental 
number characterizing the logarithmic pair (T,t). We 
have similar equations for the antiholomorphic counter- 
part (f,i. . .), with also (tt) = (Tt) = 0. 

Lattice loop models. Our aim is to measure b in per- 
colation and dilute polymers with PBC. An essential dif- 
ficulty consists in denning the models; since, for instance, 
percolation presents many non-local (or apparently non- 
local) aspects, the consideration of different geometri- 
cal features might lead to different, mutually inconsis- 
tent, field theoretic descriptions. A safe way out is to 



use supersymmetry [241 ]. but to keep things simple, we 
consider — as is quite standard — the observables encoded 
in a gas of loops, where each loop carries a weight n. 

It is well-known how to rewrite the Q-statc Potts 
model in terms of Fortuin-Kasteleyn (FK) clusters (see 



e.g. 



25|). whose hulls are dense loops with n = 
Percolation then arises as the n — > 1 limit. We choose 
PBC in both space and imaginary time directions. The 
corresponding correlations can be obtained using a trans- 
fer matrix, or, in the very anisotropic limit, a quantum 
Hamiltonian in 1+1D, H = — Y^i=i e i- where L = 2N is 
the length of the system in the space direction and the 
ei are Tempcrlcy-Lieb (TL) generators. These are repre- 
sented by the contraction e, = n of the sites i and i + 1 
on the cylinder lattice in Fig. [TJ The different sectors 
of the "Hilbcrt" space of the Hamiltonian are labeled 
by the number of through-lines 2j propagating along 
the imaginary time direction. Examples of states in the 
j = 1 sector on L = 6 sites are given by |a) = \ \ ) 
or |/3) = \y\ |^v) imposing PBC. The action of the e» 
can produce winding of through-lines around the cylin- 
der which we replace by a weight 1; note also that non- 
contractible and contractible loops have the same weight 
H [27|. 
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Under the action of H, the number of 



through-lines can only decrease, since corresponds to 
a contraction. H is thus a lower block-triangular matrix 
which corresponds in the scaling limit to Lq + Lq and we 
expect to see Jordan cells on a finite lattice. In order 
to study Jordan cells for percolation, we use the same 
prescription as in the open case (see }l4|). 

To describe dilute polymers, we use the integrable ver- 
sion of the 0(n — > 0) model on the square lattice [28| 
which yields a model of dilute loops with fugacity n = 0. 
We can write the Hamiltonian as a sum of local densities 
as before, H = — , where the operator ej is now 
more complicated (see e.g. [3]). In this case j denotes 
the number of through-lines. The states are the same as 
in the dense case, except that we now allow for empty 
sites. Examples of configurations are shown in Fig. [TJ 

Lattice Jordan cells and scalar products. Both mod- 
els have trivial partition functions, so their scaling limit 
is described by conformal theories with c = 0. The 
identity operator corresponds to the groundstate of H, 
henceforth denoted jO^), and can be found in the vac- 
uum (J = 0) sector. The state corresponding to the 
stress-energy tensor T(z) is the only state with confor- 
mal weights (h,h) = (2,0) in the vacuum sector. The 
scaling dimension A = h + h of a given excitation <j> 
with energy E^(N) can be measured numerically from 
the scaled energy gap, E^N) - E Q (N) = 2^ A + . . . , 
where Eq(N) is the energy of the groundstate and vf 
is the Fermi velocity. Bethc ansatz results 29|, 3(| yield 
vp = for percolation and vf = § for dilute polymers. 
Both the Hamiltonian and the transfer matrix commute 
with the operator u 2 that translates the sites two units 
to the right. The restriction to a given conformal spin 
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s = h — h can be done using this translational invari- 
ance. We can build eigcnspaces for u 2 and introduce a 
momentum k to characterize a given sector. It can be 
argued that the states with conformal spin s belong to 
the sector with momentum k = ^-s. Using these results, 
one can readily extract the states with conformal weights 
(h,h) = (2,0) from the spectrum. We find numerically 
that the state |TW>, which goes to |T) = T(z 0) |0) 
in the scaling limit N — > oo [up to normalization issues 
to be discussed below] , is mixed into a rank-2 Jordan cell 
with its "logarithmic partner" Note that \t^) 

belongs to the j = 2 sector. We normalize the states so 
that the Hamiltonian in the basis (|TW) , p^)) reads 



H-E (N)1 



2-kvf 
2N 



AW 2 
AW 



(4) 



Percolation 


Dilute Polymers 


L — 2N 




L — 2N 




10 


-4.33296 






12 


-4.55078 






14 


-4.68234 


10 


-4.17430 


16 


-4.76634 


12 


-4.38064 


18 


-4.82256 


14 


-4.52458 


20 


-4.86168 






22 


-4.88978 






oo 


-5.00 ± 0.01 


oo 


~ -5 


Exact 


-5 


Exact 


-5 



TABLE I: Measure of the indecomposability parameter b for 
percolation and dilute polymers. 



where A^ w 
sure b 



2N 



(T\t), we also need a 



Eq(N)). In order to mea- 
scalar product" which 



goes to the Virasoro bilinear form in the scaling limit. 
This issue was already discussed in Refs. [14, 15 1 and is 
not modified by the PBC: the natural definition of the 
scalar product for such loop models corresponds to gluing 
the mirror image of the first state on top of the second 
one. 

Virasoro algebra regularization. It could be tempt- 
ing to define a lattice version of b as the scalar product 
(t(N)\ T (N)y However, the Jordan cell in the Hamilto- 
nian is invariant under a simultaneous rescaling of the 
states \t( n >} and p w )). Therefore, we need to carefully 
normalize jT^-*) so that it goes precisely to the stress- 
energy tensor in the scaling limit 14j; since (T\T) = 0, 
the problem is non-trivial, and has in fact hindered for 
many years the numerical determination of b. Like in 
the chiral case [15j , we control normalization by using a 
lattice regularization , of the Virasoro genera 
31 1 . Introducing the notation H, 



tors 



T (N) f (N) 



arguments such as those in [3l| lead, for both our dense 
and dilute loop models, to 



N ™ 

/ 1 

ttvf 

3=1 



t inj7T /N 



12 6n >°> 



(5) 



where e M = 1 (resp. e^ = V2) for percolation (resp. 
dilute polymers) is the groundstate energy density. 

Numerical results. Gathering all these (admittedly 
elaborate) pieces, we finally give a lattice expression that 
goes to b in the scaling limit [l5| 



b (N) 



t W|#W| W 



(6) 



where we have normalized the groundstate so that 
(O^IO^) = 1. Numerical results for percolation and 
dilute polymers are gathered in Tab. [U Both theories 



seem to have b = — 5, although the precision is obviously 
much better for percolation as the Hilbert space is much 
smaller. Note that we also measured b^ in the sl(2jl)- 
supersymmetric representation of percolation [24l . [32| . 
with the same result [22j ■ Our value of b thus also applies 
to the CFT describing the Spin Quantum Hall Transi- 
tion [32 1 . 

The value b = —5 is quite surprising, since the standard 
argument identifies to a large extent the chiral sector of 
the bulk theory with the chiral theory Q, and thus one 
would expect to recover the standard values b = — | or 
|. We now turn to an analytical argument that allows 
us to derive the value 6buik = — 5 exactly. 

OPE approach. It was shown 15( that by pushing 
further the ideas in 0, Q, HH one can predict the value of 
indecomposability parameters in the chiral case using op- 
erator product expansions (OPEs). The idea is that the 
identity channel of the OPE of a primary operator $^ ^ 
(with conformal weights (h, h)) with itself is ill-defined 
as c — > 



# h , 5 (* > «)* h ,s(o > o) 



2h 



z 2 T(0) 



(7) 

where a$ is a coefficient and we focused on the holo- 
morphic sector on the right-hand side (a similar reason- 
ing applies to f). The divergence as c — >• can be 
cured if there exists a primary field X(z,z) with con- 
formal weights (hx(c),hx(c)) (2,0) that "collides" 
with the stress tensor at c = 0. One can show that a 
combination of X{z,z) and T(z) becomes the logarith- 
mic partner t(z) of T{z) at c = 0, with correlation func- 
tions like eqs. Q consistent with a value of b given by 



lim r 



c/2 



h x (c)-2 ~ nm c^o —c/2hx(c). We assumed 
that those limits exist and are equal; this turns out to be 
satisfied in the following. 

Coulomb Gas. In the remainder of this Letter, we 
shall denote the bulk fields $ r s (z) ® & r i tS >(z) using Kac 
labels, meaning that the conformal dimensions fit into 
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the Kac table as (h 
paramctrization h r s 



h r ' lS i), where we use the 



r(a:+l) — sx] 
4x(x + l) 



and c = 1 



c(x + l) 



(our models with c = thus correspond to x = 2). The 
crucial step is now to identify the field X(z,z), since 
hx(c) then determines b. In the standard approach to 
this problem 0, 0] , one looks for X among degenerate 
fields such as $3.1 or $1,5 in the chiral case. Consid- 
erations on the underlying lattice model however gives 
an unambiguous answer [33j | . First, recall that the 
dense loop model renormalizes towards a Coulomb Gas 
(CG) [34| with coupling constant g 6 [2, 4] given by 
Q = 2(l+cos — recall that n = \[Q is the loop weight. 

The associated central charge reads = 1 — 3( - g ~ 4 - ) 
and percolation corresponds to Q = 1, g = | 
larly the 0(ri) model can be expressed as dilute loops 
that rcnormalize towards a Coulomb gas with coupling 
constant g £ [1,2] given by n = —2 cos irg, so dilute 
polymers have g = §. The central charge in this case 

The operator content of the 



2ff ' 

Simi- 



reads C W = 1 - 

loop models can then be worked out from their partition 



function on a torus (24, 33 1. The candidate for the field 
X(z,z) must be a spin-2 operator, and knowledge from 
the boundary case, from our numerics, and from replica 
considerations in the case of the 0(n) model suggests 
that X(z, z) has to be a 2-leg (resp. 4-leg) operator for 
dilute polymers (resp. percolation). We find that it has 



conformal weights (h 



(2±gY 



(2-gr 



+ — ) 

T 24 I 



4g 24 ' 4g 1 24 

in both cases. Moreover this is the only possible choice 
in the CG that would yield exponents (2,0) at c = 0. In 
terms of Kac labels, we can formally identify this field as 
X(z,z) = $i._2(z) (g) $1,2(2), which is non-degenerate. 
We finally compute b using the continuation of the criti- 
cal exponents of X, finding 



lim 



c 



,{Q)/ 2 



8 (2- g )2 , c(Q)-l 



3poly 



lim 



4g 1 24 

cW/2 



I (2-g) 2 
4s 



c<")-l 
24 



-5, (8a) 
-5. (8b) 



in perfect agreement with our numerical results. 

Physical discussion. We believe that our result is 
quite crucial for the general structure of bulk LCFTs 
at c = 0. In particular, the value b = —5 should 
have profound consequences on the bulk correlation func- 
tions [3, [3^]. Indeed, the solution of the differential equa- 
tions for the four-point correlation function of the degen- 
erate bulk field $2,1(^,2:) = $2.1(2) (g> $2,1(2) (with con- 
formal weights (§, §)) or $1,3(2,2) = $1,3(2) <g> $1,3(2) 
(with conformal weights (i, i)), and the resulting confor- 
mal blocks lead unambiguously to the values 6 po i y = f 
and 6p 0r cn = — | respectively (this is one of the arguments 
used in |4| to predict the value of b for percolation and 
polymers) , provided that one ignores normalization issues 
in the bulk. That b takes none of these values forces one 



to conclude that, within the c = conformal field theory 
describing percolation (resp. polymers), the four point 
correlations of the energy operator $2,1 (resp. $1,3) 
must be set equal to zero exactly. This is compatible 
with the known — but maybe underestimated — fact that 
one cannot normalize these correlation functions to ren- 
der them non-trivial and without divergences [H. Note 
also that we were able to show from both algebraic and 
numerical arguments that the energy operator $2,1 (resp. 
$1,3) is mixed into a rank-2 Jordan cell, and hence has 
a zero norm-squared. One can deduce that all its cor- 
relation functions are actually vanishing |37| . and thus 
irrelevant for the determination of b. Note that the van- 
ishing of the multi-point correlation functions of the en- 
ergy does not preclude other correlators involving the 
energy from taking non-zero values. The computation 
of non-vanishing correlation functions in general never- 
theless appears unreachable. Indeed, since b = —5 and 
occurs, from our OPE approach, because of mixing with 
operators having negative Kac labels, progress would pre- 
sumably require us to control the four-point functions 
of such operators — a situation for which the differential 
equation or CG approaches are known to fail. 

Note that the crucial point of our calculation is the 
proper identification of the field X(z, z) among primary 
operators given by the CG. We emphasize that within 
this natural "minimal" geometrical formulation of the 
problem, there is no bulk counterpart of the boundary 
fields $3,1 or $1,5 that played the role of X in the chi- 
ral case. Instead, we have shown that X(z, z) must be a 
non-Kac operator with negative Kac labels, thus yielding 
a value of b different from what was expected from, e.g., 
Refs. [1,[35{. The possibility of having a value of b in the 
bulk different from the boundary case was suggested in 
Ref. [35j]. We emphasize that our geometrical bulk the- 
ory admits a unique value of b [40T ] . Note however that a 
"non- minimal" formulation — e.g., involving further geo- 
metrical observables [HI, |39| — might eventually lead to a 
more complicated Jordan structure for T(z). Finally, we 

remark that -r^ — = 7-^ h , 1 ; this is an intriguing ob- 

servation indeed, since b~ x is the parameter that appears 
multiplicativcly in correlation functions and OPEs. 

In conclusion, we see that the complete picture of the 
bulk theory of percolation (or polymers) is more compli- 
cated than expected, cannot be seen as a simple gluing of 
two chiral theories, and seems in fact to have very little 
to do with the nice understanding gained over the last 
few years in the boundary (chiral) case. Further under- 
standing follows from a detailed algebraic analysis, which 
will appear elsewhere (22j . 
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